Nonmonotonous magnetic field dependence and scaling of the 
thermal conductivity for superconductors with nodes of the order 

parameter 

Yu. S. Barash a ' and A. A. Svidzinsky a,c 
a) P.N. Lebedev Physical Institute, Leninsky Prospect 53, Moscow 117924, Russia 
b) DRFMC, Service de Physique Statistique, Magnetisme et Supraconductivite, CENG, 
17, rue des Martyrs, 38054 GRENOBLE - Cedex 9, France 
c) Department of Physics, Stanford University, Stanford CA 94305-4060, USA 

(February 1, 2008) 

Abstract 

We show that there is a new mechanism for nonmonotonous behavior of 

magnetic field dependence of the electronic thermal conductivity k of clean 

superconductors with nodes of the order parameter on the Fermi surface. In 

particular, for unitary scatterers the nonmonotony of relaxation time takes 

place. Contribution from the intervortex space turns out to be essential for 

this effect even at low temperatures. Our results are in a qualitative agreement 

with recent experimental data for superconducting UPt% MM. For i?2u-type of 

pairing we find approximately the scaling of the thermal conductivity in clean 

limit with a single parameter x = ^-y at low fields and low temperatures, 

as well as weak low-temperature dependence of the anisotropy ratio K zz /K yy 

in zero field. For £i ff -type of pairing deviations from the scaling are more 

noticeable and the anisotropy ratio is essentially temperature dependent. 
PACS numbers: 74.25.Fy, 74.70.Tx 
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I. INTRODUCTION 



Growing amount of experimental data indicates that many of high-temperature and 
heavy-fermion superconductors have highly anisotropic order parameter on the Fermi sur- 
face. The possibility for the order parameter to have opposite signs in different regions on 
the Fermi surface, lines of nodes between those regions (and, possibly, point nodes too) 
has attracted much attention. Since the presence of nodes of the order parameter leads to 
the existence of low-energy excitations in spectrum, this strongly modifies low-temperature 
behavior of thermodynamic and transport characteristics both in the absence and under the 
applied magnetic field ( see reviews Hf| and for more recent literature, for example, |5|-|8| ). 
For these reasons measurements of temperature and magnetic field dependences of those 
characteristics are important experimental tools for probe the anisotropic structure of the 
order parameter 12j4^|. As compared to thermodynamic characteristics like specific 

heat, an important advantage of studying transport properties, in particular, the thermal 
conductivity is that it is a directional probe, sensitive to relative orientations of the thermal 
flow, the magnetic field and directions to nodes of the order parameter. 

For isotropic s-wave superconductors low temperatures are obviously defined as satisfying 
the condition T <C A max , when the number of quasiparticles thermally activated above the 
gap is exponentially small. In the presence of nodes of the order parameter on the Fermi 
surface an additional small energy scale 7 appears, describing the bandwidth of impurity- 
induced quasiparticle bound states. Then the low-temperature region should be divided 
into two parts. For temperatures T < 7 transport properties are dominated by bound 
states. Non-zero density of these states on the Fermi surface results in linear temperature 
dependence of the thermal conductivity in this case. At the same time under the condition 
7 ~ T <C A max , which can be satisfied for sufficiently clean superconductors, one can 
disregard the influence of bound states. Then the presence of nodes of the order parameter 
leads to characteristic higher order power-law behaviors of the thermal conductivity with 



temperature [H-21 
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The thermal conductivity of s-wave isotropic superconductors in the presence of the ap- 



plied magnetic field exhibits nonmonotonous dependence upon the field p2|-f25|j. In high 
fields the thermal conductivity k(B) rises rapidly due to suppression of the order parameter 
as the magnetic field approaches At low fields k(B) decreases with increasing mag- 
netic field. This phenomenon is usually attributed to the scattering of electrons (or/and 
phonons) by the fluxoids, which was theoretically considered in [pq| . However, in isotropic 
s-wave superconductors at low temperatures the number of quasiparticles thermally acti- 
vated into scattering states are exponentially small, so that the contribution to the thermal 
conduction along the magnetic field from bound states within the vortex cores become es- 
sential. In contrast to isotropic superconductors, for anisotropically paired superconductors 
with nodes of the order parameter on the Fermi surface, the intervortex space can domi- 
nate thermodynamic and transport characteristics at low temperatures due to low energy 



excitations with momentum directions near the nodes ||[7|,^7|-|29| . Under these conditions 



the characteristics usually become quite sensitive to the applied magnetic field even if it is 
of relatively small value. For example, quasiparticle density of states on the Fermi surface 
for clean superconductors with nodes, takes nonzero value just due to the applied magnetic 
field ( as well as impurities ). Since low-temperature behavior of thermodynamic quantities 
like the specific heat is directly associated with the behavior of the density of states on the 
Fermi surface N s (0), this is the reason for their sensitivity to the magnetic field |5]||. For 
transport characteristics additional quantities are of importance in this respect. These are 
scattering relaxation times for various channels. For instance, it is known that the relaxation 
time t s (u) at sufficiently low energies for elastic impurity scattering in clean superconductors 
with nodes, may be quite small for unitary scatterers and extremely large in the case of Born 
scatterers 0,^1|,|T^,^] . While t s (u) for Born scatterers is directly associated with N s (cu), 



in the unitary limit t s (uj), as a function of energy, is not reduced entirely to the density of 
states and should be considered as an independent quantity. 

Below we show that for anisotropically paired clean superconductors with nodes of the 
order parameter on the Fermi surface, at low temperatures there is important additional 



mechanism for nonmonotonous dependence of the electronic thermal conductivity upon the 
magnetic field. This mechanism is associated with the electronic contribution to k(B) mostly 
from the intervortex space due to the presence there of extended quasiparticle states of 
low energies with momentum directions near nodes of the order parameter. Contrary to 
isotropic s-wave superconductors, the main effect comes in this case from the influence of 
condensate flow field (even of relatively small value) on the quasiparticle impurity scattering 
in the unitary limit, rather than from scattering of quasiparticles directly on vortex cores. 
For these two types of scattering inverse relaxation times could be added under certain 
conditions ||26|| . However, the contribution from scattering of quasiparticles by vortex cores 
is supposed to be negligibly small under the conditions considered below, as compared to 
the one from scattering by impurities in the intervortex space. 

In the unitary limit relaxation time for scattering by nonmagnetic impurities of low- 
energy quasiparticles, is found below to be a nonmonotonous function upon the condensate 
flow field. This takes place even for uniform superfluid flow, that is in the absence of 
any scattering by vortices. Furthermore, for type II superconductors with large Ginzburg- 
Landau parameter, the superflow induced by magnetic field in the intervortex space can be 
considered on sufficiently large distances from vortex cores (much greater than the coherence 
length) as quasihomogeneous flow. This allows one to consider approximately the thermal 
conductivity as a function of the local value of condensate flow field, as this would be for the 
uniform flow. Magnetic field dependence of the thermal conductivity obtained in the first 
approximation as spatial averaging of the result over the intervortex space, is found below 
to be nonmonotonous under certain conditions. 

This effect can be important, in particular, for the analysis of recently observed non- 
monotonous behavior of the thermal conductivity in UPt 3 , since the electronic contribution 
to n is known to dominate there below IK . We consider (1, i) phases both for the E 2u 
representation and for the E\ g one as candidates for the type of superconducting pairing 
in UPt 3 at low temperatures and under the weak applied magnetic field B c \ < B B c2 
(see, for example, |3|J^j3^|). Our theoretical results, basing on this consideration, are in a 
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qualitative agreement with those experimentally observed in Mostly the both models 

give rise to such an agreement. The difference between corresponding predictions is not too 
great, although it is of importance permitting to distinguish between them. Recent exper- 
imental data i n particular, allow for determining the low-temperature interval, where 
the power law temperature dependence of the thermal conductivity takes place. We find, 
that the behavior of the zero-field anisotropy ratio for the thermal conductivity for these 
temperatures seems to indicate in favor of E2 U type of pairing. At higher temperatures 
the behavior of the thermal conductivity becomes essentially depending upon the particular 
form of the order parameter all over the Fermi surface, not only near its nodes. Under this 
condition there are various possibilities to fit experimental data within the framework of 



both models, so that the problem to distinguish between them becomes ambiguous one |]T8 
By contrast, at sufficiently low temperatures T T c the behavior of the thermal conduc- 
tivity is governed mainly by the behavior of the order parameter near nodes, as well as by 
the strength of scatterers. This leads, in principle, to the possibility to identify the behavior 
near the nodes and, hence, the type of superconducting pairing, but not a particular form of 
the order parameter all over the Fermi surface. This circumstance was already emphasized 
earlier in [HJ, where the accent was made on the ultra low temperatures T < 7. Quite a 



small value 7 ~ 0.017/C taking place for clean samples of Ref. give the possibility to 
determine and concentrate on the low-temperature region 7 ^ T < T c , while the ultra low 
temperatures in this case seem not to be sufficiently studied experimentally yet. 

The article is organized as follows. Basic equations are listed in the next section. 
We consider clean superconductors at low temperatures under the applied magnetic field 
( H c \ ^5 B <C H c2 ) and disregard the contribution from impurity bound states up to Sec. 6. 
In Sec. 3 we determine the nonmonotony of the relaxation time in the unitary limit. It 
is shown, that for the energy dominating regime PfV s <C uj the relaxation time t s (uj,v s ) 
decreases with increasing superflow in the presence of the line of nodes of the order pa- 
rameter, in contrast to the case of superflow dominating condition u <C PfV s , when the 
relaxation time rises. Analytical calculations are carried out for sufficiently small quantities 



PfV s , oj <C A max . In Sec. 4 the thermal conduction as a function of the superflow as well 
as along the applied magnetic field is analyzed. We examine three important examples of 
anisotropic gap functions: the polar state and (l,i) superconducting states for E\ g and E 2u 
types of superconducting pairing in hexagonal crystals ( like UPt^ ) . We find that in the 



temperature dominating region T C J~§— < T < A m(M . the thermal conductivity diminishes 



with the magnetic field for all three types of pairing considered, if the magnetic field is 
applied within the basal plane. This turns out to be valid as well for k zz and the magnetic 
field along z-axis (of high symmetry) for E\ g and E 2u types of pairing, not for the polar 
state. The contribution to the magnetic field dependence of k from the intervortex space 
r ^> £ dominates in this case, which justifies the semiclassical approximation we use. For 
higher fields or lower temperatures the thermal conductivity always rises with the increasing 
magnetic field, although in this case intervortex space dominates only for the magnetic field 
oriented within the basal plane. This ensures the nonmonotonous magnetic field dependence 
of k ( at least for E 2u and Ei g models ) under the condition @T C J < T <C T c , where 
numerical factor f3 may be of the order of unity. In Sec. 5 the low-temperature behavior 
of the thermal conductivity in zero field and the magnetic field dependence of k,u(T,B) at 
low fields and low temperatures are studied in clean limit (the z-th axis is aligned along the 
magnetic field). We show for the temperature region 7 < T <C A max , that scaling of the 
thermal conductivity with a single parameter x = if-\] & t low fields and low temperatures 
(both for i = y and for % = z) as well as weak low-temperature dependence of the anisotropy 
ratio K zz /Kyy in zero field, are approximately valid in clean limit for (l,z)-phase of E 2u -type 
of pairing. Qualitatively it is quite close to what it is observed experimentally for UPt^ in 
]T|Q]. Under the same conditions E\ g model results in more noticeable deviations from the 
scaling, and in essential temperature dependence of the ratio k zz (T,B = 0)/ K yy (T, B = 0). 
New test is suggested for discrimination between candidates for the type of pairing in UPt^, 
based on the dependence of k zz upon the value of transport supercurrent flowing in thin 
films or whyskers along the hexagonal axis in the absence of the magnetic field. In Sec. 6 
we describe the behavior of the thermal conductivity under the condition PfV s , T < 7, when 
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contribution from impurity bound states dominates the thermal conductivity. We find in 
this limit, that the thermal conductivity monotonously decreases with increasing condensate 
flow field (at least for 7 <C A max ) and does not satisfy the above-mentioned scaling behavior. 
This is also in agreement with the experimental results for UPt$. 



II. BASIC EQUATIONS 

Taking into account influence of the magnetic field (satisfying the condition B c \ < B <C 
B c2 ) on the thermal conductivity, we are interested in the contribution from large distances 
from vortex cores (£(T) < r), where the problem can be considered approximately as locally 
quasihomogeneous one on the basis of semiclassical approximation ||. Being justified for 
those cases, when the contribution from the intervortex space £ ^5 r turns out to be domi- 
nating, such an approach simplifies greatly all analytical considerations permitting to obtain 
correct results up to numerical coefficients of the order of unity. Besides, we assume the 
particle- hole symmetry and consider superconducting states to be unitary: AA^ = | A| 2 cr , 
where <3"o is the 2x2 unit matrix. Then the expression for the thermal conductivity under 
the applied magnetic field can be written in the form 

13 2 Jo \TJ cn , h 2 f^_\lmu\ 2{Re Co -v f v s ) 



277 



\Cd - v f v s \ 2 - {Cb- VfV s ) 2 



xRe 2 v /(^-^v g r-|A(p / )P+ l ^ S| ; - ' \) ■ (!) 

Here v s = (1/2) (V% + (2e/c)A) is the gauge-invariant condensate flow field induced by the 
magnetic field (see, for example, x is the phase of the order parameter) , Co and A p 



are the quasiparticle energy and the order parameter renormalized by impurities; notation 
(. . .) Sf means the averaging over the Fermi surface. Further we ignore contributions to the 
quasiparticle group velocity due to the momentum direction dependence of the order param- 
eter, since for small enough parameter A max /ef this would lead to a quite small corrections 
to the thermal conductivity ||34j| . Also we assume below for simplicity a spherical Fermi 



surface and introduce the superfTuid velocity v s = v s /m, so that under this simplification 

VfV s = mVfV s = PfV s . The branch of square root function in (|1|) should be chosen to have 

nonnegative imaginary part. 

Considering sufficiently clean superconductors in the temperature region 7 < T, one can 

disregard the contribution from impurity bound states. Then Eq.(|IJ) reduces to 

«ij 9 fu\ 2 1 t s (lu,v s ) 

duJ \^) —TTT\ Z %(^>0 • ( 2 ) 



k n {T c ) 2vr 2 T c io \T) ^2 (^_\ r N 

Here kn(T c ) is the thermal conductivity in the normal state at T = T c ; r a (u) and are 
relaxation times for quasiparticle scattering on impurities in superconducting and normal- 
metal states respectively. The quantity 1^(0), v 8 ) in the case of spherical Fermi surface is 
defined as follows 

|A(p/)|<|w-p/v«| 



I i:j (uj,v s ) 



,,1- ' A ^ 2 . (3) 



Assuming particle-hole symmetry, one can represent the relaxation time in the Born 
approximation in the form: 

T a (u) 1 



r N ImG (w) 
while in the unitarity limit the corresponding expression is 

T s (u) _ 1 

r N Im G (u) 



(4) 



(5) 



The quantities G (u>), G\(uS) are 

n ( \ ■ f dQ u ~ Pf Vs rets 

G oH = % \ — 1 , (6) 

Gi(u) = i — v 7 ; (7) 

J ^ ^J(uj-p f v s ) 2 -\A(p f )\ 2 

In the absence of the superfluid velocity above expressions are reduced to well-known rela- 
tions |TJ]]. For the chosen branch of square root function the real part of the integrand in 
(|) (that is the density of states) is a nonnegative quantity. 
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III. RELAXATION TIME IN THE UNITARITY LIMIT 
A. Magnetic field, aligned along z-axis 

Let us consider superconducting phases with order parameters having line of nodes in 
the equatorial plane on a (spherical) Fermi surface, and the applied magnetic field satisfying 
the condition B c \ < B B c2 and aligned along the crystalline axis of high symmetry 
(z-axis). As particular examples, there may be the polar phase and (1, i)-superconducting 
phases of E\ g - and i^u-types of pairing in hexagonal superconductors. Last two phases 
are widely discussed as candidates for the type of superconducting pairing taking place 
in the heavy fermion superconductor UPts P,|^J5^]. All these particular order parameters 
change their signs under the reflection across the equatorial plane on the Fermi surface: 
A(tt — 9) = —A(9). Due to this property function Gi(u) is equal to zero for the magnetic 
field parallel to z-axis. 

The applied magnetic field influences the relaxation time mainly due to its dependence 
upon the superflow, induced by the field and aligned in the particular case parallel to the 
basal plane. In calculating the dependence of the relaxation time t s (uj) upon the superflow 
under the condition PfV s , u <C A max , it is essential that only narrow regions near nodes of the 
order parameter on the Fermi surface govern the dependence of Gq(uj, v s ) — ReGo(u, v s = 0) 
upon the superfluid velocity. At the same time the quantity ReGo(u>,v s = 0), generally 
speaking, is formed by the whole Fermi surface. We assume that the order parameter near 
the line of nodes (|~ — 9\ <C 1) may be represented as |A| = AjJ^ — 9\, while near the point 
node on one of the poles (\9\ <C 1) it has the form |A| = A2|#| r \ n = 1,2. Calculations show, 
that if one keeps terms up to p/f s /Ai j2 in the expansion in powers of the superfluid velocity, 
taking into account as well the first power of small parameter u;/A li2 , then contributions 
from the line of nodes and quadratic points are essential. At the same time contributions 
from linear points may be neglected as compared to ones from the line of nodes. As a result 
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we obtain the following expressions for the relaxation time in the unitarity limit 



r,U') ( ^ ^ Ai 



/ 

+ ln 



2Ac 



2cj 



vrAx 1 + 



2A 9 



Ai 



ReG (o;,u s = 0)+ 



2uj 



ye Iuj + ^uj 2 - (p f v s ) 2 j J 



UJ 



PfV s < uj < A 

7: 



(8) 



r a (a?) _ p f v s 



Ax 



\ 



UJ 



+ 



Sill 



K PfV s j p f V s 



K PfV s 



tiuj 
4AT 



uj 



2 * Al ReG (cJ,^ = 0) + ln' 2w 



+ 



UJ 



K PfV s 



+ 



sin 



7ra;Ai 
4:p f v s A 2 



uj < Pjv s A r 



(9) 



The most important feature of the dependence upon the superflow, describing by Eqs.@, 
(§), is its nonmonotonous behavior. Indeed, one can easily see under the conditions 
p/v s <C uj <C A max , that the relaxation time decreases with the increasing superfluid ve- 
locity according to the relation 

2 



T S (UJ,V S ) 7TUJ ( i Ai 



2Ai 



2A, 



7T 1 + 



'A 

. , , — (ReG (^,^ s = 0)) 2 - 
^i \ V uj 



2A, 



-^0^ReG o (uj,v s = O)\ . (10) 

At the same time under the opposite condition uj <C pfV s <C A max we see from @, that 
r s/ r iv = P/^s/Ai increases linearly with the parameter p/u s . 

Expressions (|8|)- ([T0|) describe directly the relaxation time for superconducting state 
of i?2u-type of pairing, for which the order parameter has the line of nodes on the equator 
and quadratic point nodes on the poles. Results for Ei ff -type of pairing ( line of nodes on 
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the equator and linear point nodes on the poles ) and for the polar phase ( nodes are only 
on the line on the equator ) follow from fl5|)-(|lCJD in the limit A 2 — > oo. 

In calculating KeGo(uj,v s = 0) one should describe the particular form of the order 
parameter all over the Fermi surface. In the particular case of polar phase A = Aicos# 
explicit integration in (|6[) results in the following expression 



ReGnf^, v* 



UJ 



Ai 



A: 



UJ 



Ax 



UJ 



(11) 



Substituting (|TTD into (||) and (Q), one gets that at the point of minimum of the relaxation 
time the superffuid velocity satisfies the condition PfV s ~ uj\n (Jio^- Under the condition 
PfV s <wC Ai relaxation time decreases with increasing parameter pfV s : 



t s {uj,v s ) 



TN 



2uj 
ttAi 



7T~ 



+ In' 



/2Ai 



V uj 



2a/ 



(12) 



while in the opposite case w < p/u s < A] we have 



T S (UJ,V S ) _ p f V s 



TN 



Ax 



0J 



In 2 



' 4Ai 

ep f v s 



(13) 



If coin (^j) ~ Pf v s, expression (p~3|) increases with increasing pfV s . 

One can see, that nonmonotonous behavior of the relaxation time with superfluid velocity 
takes place in the unitary limit due to different behaviors of ReGo and ImGV At sufficiently 
low energy the density of states (and, hence, ImG ) increases with increasing magnetic field, 
while ReGo turns out to decrease. In the case in question the relaxation time is described 



by the relation [35 



^ = ImG„( W ) + (ReGo(w)r 



r N lmG (uj) 
and may manifest nonmonotonous magnetic field dependence (see Fig. 1). 



(14) 



B. Magnetic field parallel to the basal plane 



Nonmonotony of the relaxation time as a function of the applied magnetic field may take 
place for various orientations of the field, due to the same reasons as for the field parallel to 
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the z-axis. In this subsection we present shortly some results on magnetic field dependence 
of relaxation time for the field lying within the basal plane for polar phase (A(pf) = A p (6)), 
(l,i)-phase of £'i g -representation of D Q h point group (A(p/) = Ai g {6) exp(iyj)), and (l,i)- 
phase of ^^-representation (A(p/) = A 2u (6 l ) exp(2i</?)). Spherical angles 6 and y> character- 
ize here momentum directions on a (spherical) Fermi surface, and functions A p (9), Ai g (8), 
A2u(Q) are supposed to be real. 

Let ?/-axis be parallel to the magnetic field. We consider some fixed spatial point far 
enough from the vortex core, where induced superfluid velocity, lying in xz-pl&ne, constitutes 
angle 9q with z-axis. Although for the field orientation within the basal plane Gi(u>) differs 
from zero, the analysis shows that this quantity may be neglected as compared to Gq(uj), 
at least for the types of pairing considered and under the condition u, pfV s <C A max . This 
allows one to use Eq.fll4]) in calculating t s (uj). Assuming the same relations for the order 
parameter near the line of nodes (|A| = Ai|| — 9\) and point node (|A| = A 2 |0| n , n — 1,2) 
as earlier, we get again (within the same accuracy) that the relaxation time diminishes with 
increasing superfluid velocity under the conditions pfV s <C uj <C A max : 



T S (U,V S ) _ TTUJ f 1 _Aj_ 



2Ai 



2A 5 



7T 1 



2A-, 



— (ReG (uj,v s = 0)y 



UJ 



1 + a; 



A A (j)fV s ) 2 cos 2 6> c 



2a/ 



ReG (^, v s = 0) 



In the opposite case uj -C pfV s <C A max relaxation time increases with v s : 



r s (v s ) 



PfV s 



sin | tt \ cos 6 \ 



(15) 



(16) 



r N - \ Aj 4A 2 J 

Eqs.(|l5D, C p76|) describe the relaxation time for (l,z) superconducting state of E^u-type 
of pairing. As it was in the case of Eqs.(|8|)- ([10|) , results for E\ g -type of pairing and for the 
polar phase follow from (|l^), (|i~6"D in the limit A 2 — > oo. For the polar phase we obtain 
under the condition uj <C A\ the explicit expression ReG (u,v s = 0) = (uj/ Ai) ln(2Ai/a;). 
Then Eq.(jl5|) reduces to 



r s {u,v s ) 



2uj 

7lA 1 



4 \ UJ 



1 ( pfV s cos 9 \ / 2Ai 

2 I uj ) n V UJ 



(17) 
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So, magnetic field dependence of the relaxation time for the magnetic field oriented 
within the basal plane, turns out to be qualitatively quite close to the one for magnetic field 
orientation along the z-axis. 



IV. THERMAL CONDUCTIVITY 

In this section we consider the electronic thermal conduction along the direction of the 
applied magnetic field. Let the magnetic field be aligned firstly along the z-axis, so that the 
supercurrent flows parallel to the basal plane. 

For low temperatures and small superfluid velocities satisfying the conditions T, p/v s <C 
A max , there are two characteristic regions pfV s < T <C A max and T <C pfV s <C A max , 
where dependences of the thermal conductivity upon v s may differ essentially. While within 
the region T <C PfV s <C A max the thermal conductivity always rises with increasing v s , it 
may be possible both increase or decrease under the condition pfV s < T < A max . Thus, 
nonmonotonous behavior of the relaxation time may result in respective nonmonotonous de- 
pendence of the thermal conductivity upon the superflow (and eventually upon the magnetic 
field) . However, nonmonotony is not a inevitable consequence for the thermal conductivity, 
which would be valid for any types of pairing in the presence of such a behavior of the 
relaxation time. The point is that the quantity Iij(u,v s ), which forms the behavior of k 
along with the relaxation time (see (0)), usually turns out to be monotonously increasing 
function of v s . As a result, the product t s (cj, v s )Iij(u, v s ), containing nonmonotonous func- 
tion t s (u,v s ), may manifest both monotonous or nonmonotonous behavior, depending on 
the particular behaviors of r s (u,v s ) and Iy(uj,v s ). 

For sufficiently low temperatures T <C p/v s <C A max (or, in other words, for large enough 
superfluid flow), it follows from (0) in the first approximation 

Kij 3Tt s (uj = 0,v s ) 

TtF\ = t Z h 3 (uj = 0,v s ) . (18) 

k>n{J-c) J-c Tn 

The possibility to disregard the contribution from impurity induced bound states in this 
case is justified under the condition pjv s ^> 7. 
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One can show, that the contribution to I zz (uj = 0, v s ) from narrow region around the line 
of nodes is (pfV s ) 3 / 12 Af. Linear point nodes do not contribute to I zz (u = 0, v s ), so that 
Eq.(|18]) is not the appropriate approximation for the case of the order parameter having only 
linear point nodes. It is important, however, that the contribution from linear point nodes 
in this limit is evidently less than from the line of nodes. At the same time the contribution 
from the quadratic point nodes dominates being equal to (pfV s ) 2 /QA 2 .. Furthermore, since 
ReGo(^ = 0,v s = 0) = 0, we find from @: t s (u = 0,v s )/tn = PfV s / Ai. The contribution 
to t s (uj = 0,v s )/tn from quadratic point node is (pfV s ) 2 /2A|, which can be disregarded as 
compared to the above value of t s (uj = 0, v s )/tn formed by the line of nodes. Subsequently, 
under the condition T <C pfV s <C A max we find in the case of polar phase and for (l,i)-phase 
of Eig-type of paring: 

T (19) 



k n (T c ) 4T C V A 1 
For (l,i)-phase of E 2u -tjpe of paring it follows, that 

= JL^sf 

k n {T c ) 2T C A x Aj y 1 

These higher power law dependences upon the local superfluid velocity result after av- 
eraging over the intervortex space in linear dependence upon the magnetic field, formed 
by the lower border r ~ £ of the intervortex space considered. Indeed, superfluid velocity 
v s [r) depends upon the distance from the vortex core and in the simplest case of circu- 
lar supercurrents one has v s (r) = l/2m e r for ( C r < A. Then we obtain after spatial 
averaging: 

Ku(T , H) ^(^y^ Ks{Tc) (21) 

for polar phase and for (l,i)-phase of E^-type of paring, and 

MT.^k^^MT.) (22) 
for (l,i)-phase of -E 2tl -type of paring. 
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Since regions close to vortex cores (r ~ £) dominate in spatial averaging of the thermal 
conductivity under the particular conditions T <C p/v s , H\\Oz, the approximation of locally 
quasiuniform superflow, generally speaking, is not directly applied to this case. In addition, 
the relation PfV s (£) ~ T c guarantees validity of inequality T <C pfV s for all T < T c and, 
from the other side, violates the condition pfV s <C A max , which leads to the conclusion 
about rough qualitative character of the particular results derived above. Strictly speaking, 
in this particular case it is necessary to check the contribution to k not only from the 
delocalized quasiparticles, but from quasiparticles bound to vortex cores as well. That 
problem is, however, beyond the present article. It is important, that in any case the 
thermal conductivity increases with the increasing magnetic field in the limit considered. 
This is, fortunately, the only particular example in the article, where the intervortex space 
doesn't dominate. Furthermore, it is known for s-wave superconductors, that contribution 
from bound states ( localized within vortex cores ) to the low-temperature behavior of 
the thermal conductivity along the magnetic field, is negligibly small as compared to the 



normal state thermal conductivity [23J. In the presence of nodes of the order parameter 
extended (scattering) low-energy quasiparticle states contribute much greater. Indeed, as 
p f v s (£,)/A max ~ 1, we get from Eqs.(|2T|), (f^) k zz {T) ~ {B / B c2 )k n (T). This estimate is 
quite close to that observed experimentally for UPt^ for low temperatures fl],|2| and about 
on two orders greater, than the one realized for conventional superconductors ||24)| . Also we 



note, that for the field within the basal plane large distances from the vortex cores dominate 
for any relation between T and PfV s (see below) ensuring the applicability of the semi classical 
approximation in that case both in temperature dominating regime T 3> PfV s and under the 
superflow dominating condition T -C pfV s for three types of pairing discussed. At last, for 
Born scatterers the situation changes and in spatial averaging in the superflow dominating 
limit large distances from vortex cores dominate both for k zz and the field along the z-axis 
and for K yy and the field along the y-axis. 

Let us discuss now the other case pjv s « T « A ma2 .. It is worth noting that this 
condition can't be valid close to vortex cores (for r ~ £). However, as we show below, under 
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the condition p/v s < T < A maa; fairly large distances from vortex cores (of the order of 
the intervortex distance R ~ £y^^ ^ ma y dominate in the spatially averaged thermal 
conductivity. 

Considering firstly the thermal conductivity for the polar phase and substituting A = 
Ai cos 9 into (|3]), we obtain under the condition u, pfV s <C Ai : 



LO 



[2u 2 + 3(p f v s ) 2 } sin" 1 



LO 

p f v t 



+ [3^ 2 + 



16A? 



+2(p f v s ) 2 ] sJ{pfV s ) 2 -u 2 - | \{p f v s ) 2 - u 2 f /2 , lo < p f v s 



(23) 



7TL0 3 + ^-Uj(p f V s ) 2 



LO > PfV s . 
'PfV s 



For to < pfV s it follows from (H), ©, that I(u,v 8 ) < [^f-J , r s ~ HJr)- Taking 
account of these estimates, we find, that under the condition PfV s < T < Ai the contribu- 
tion to the thermal conductivity (0) from the integration over the region to < pfV s turns out 
to be ^ jffi^/^ • Since we are interested in contributions to k up to (pfV s ) 2 , one should 



consider the integration over to > pfV s and substitute into @ equations (|12|), (|23|): 

k« 161. 74T 5 
^(T c ) ~ T c Aj 



In 2 r^-1.631n^ 1+3.16 
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2 ^2Ai 



- 1.67 hi [ ''hjT' ) 3.1s 



(24) 



We see from (pTj) and (0), that k 22 for polar state is a monotonously increasing function 
of the superfluid velocity. According to ( p4[) the thermal conductivity for polar state in the 
presence of magnetic field roughly oc T 3 , while in the absence of magnetic field oc T 5 . 

By contrast, the thermal conductivity of (l,z)-phases for E\ g and E 2u types of pairing 
turns out to be nonmonotonous one. The point is that the contribution to I zz from linear 
point nodes is approximately to 2 /3A 2 ,, while from quadratic point nodes is 7tlo/8A 2 . These 
terms dominate as compared to the contribution from line of nodes. This means that 
magnetic field dependence of the thermal conductivity is formed in the cases in question 
entirely by the relaxation time: 
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(25) 



This nonmonotony is seen in Fig. 2 where we plot the thermal conductivity for three 
particular types of pairing and various values of temperature as a function of the superffuid 
velocity |59fl . 

After averaging over the intervortex space, where expressions (^4j), (|2~5"D are valid, we 
obtain 



Sign plus corresponds here to the case of polar phase, while sign minus - to E\ g and E 2u 
types of pairing. Since ReGoiTu, 0) is roughly proportional to T (neglecting the logarithmic 
factors; see below Eq.(|33D), positive dimensionless function A ZZ (T) turns out to have a linear 
temperature dependence for E 2u model, a quadratic temperature dependence for Ei g , while 
for the polar phase A ZZ (T) is roughly a cubic function of temperature (see Eq.(|2~3|) as well). 
Large logarithmic factor in magnetic field dependence (|2B|) is formed by contributions from 
sufficiently large distances from vortex cores (of the order of intervortex distance R). This 
justifies the approach we use, which is based on the consideration of quasihomogeneous 
superconducting intervortex space and valid for large distances from vortex cores r ^> £. 

So, in spatial averaging of the thermal conductivity the intervortex space, generally 
speaking, is divided into two parts. The first one is defined by the condition T < PfV s {r) 
and contributes to the monotonously increasing thermal conductivity with the magnetic 
field. The other part of the intervortex space pfV s (r) < T includes fairly large distances 
from vortex cores. For Ei g and E 2u representations it contributes to decreasing thermal 
conductivity with increasing magnetic field and turns out to be essential (that is may com- 
pete with the contribution from the first part) only under the condition pfV s (R) -C T. This 
condition reduces to the inequality T ^> T C ^J -J^ which has to be satisfied ensuring the pos- 
sibility for the nonmonotonous magnetic field dependence of the thermal conductivity. Thus, 




(26) 
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in the presence of nonmonotony the minimum of k(B) has to be shifted gradually to lower 
fields with decreasing temperature. Below some characteristic temperature ~ /3T c y^l this 
nonmonotony disappears at all. Rough estimates show that numerical factor (3 may be of 
the order of unity. This is in agreement with that observed experimentally in UPts [pl0. 



Naturally, the relation (3T c ^^j^- < T < T c can be satisfied only for large enough value of 
the Ginzburg-Landau parameter as this takes place, in particular, in UPt^. 

Let now the applied magnetic field lie in the basal plane along y-axis. While the re- 
laxation time behavior is quite similar for both orientations of the magnetic field, it is not 
always the case for Jy. As a consequence, a qualitative behavior of the thermal conduction 
along the magnetic field for B\\Oz and B\\Oy may differ as well. In particular, we show 
below for the field orientation B || Oy, that nonmonotonous magnetic field dependence of 
the thermal conductivity K yy takes place for all three types of pairing in question, including 
the polar phase (in contrast to k zz for the orientation B||Oz). 

Under conditions oj, pjv s <C A max we get, that the main contribution to I yy comes only 
from the narrow region of the line of nodes for all types of pairing considered (including 
E 2u - representation with quadratic point nodes on the poles of the Fermi surface): 
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(27) 



PjV s sm & V 

As it follows from (P7[), l yy {oj = 0,v s ) = pfV s \sm9 \/6A 1 . Substituting this value 
together with into flT£|), we obtain the thermal conductivity at low temperatures 
T pjv s A max for (1, z)-phase of E 2u -tjpe of pairing: 
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(28) 



As earlier, the answer for polar phase and for (1, i)-phase of _Ei 5 -type of pairing is obtained 
from (|28|) in the limit A2 — > 00. 
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For fairly low superfluid flow, when the condition pjv s < T <C A moa; is satisfied, the 
expression for the thermal conductivity can be obtained, substituting I yy = -g^- together 
with Eq.pl) into Eq.(0): 
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In the particular case of polar phase we obtain from 



11) the following dependence 



of the thermal conductivity upon the low enough superfluid flow: 
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(30) 



Spatial averaging of K yy over that part of the intervortex space, where Eq.(^) is valid 
( that is over its exterior part pfV s (r) < T ), leads to the thermal conductivity, which 
diminishes with increasing weak magnetic field for all three types of pairing according to 
Eq.(p6|) (where one should change indices z — > y and retain only sign minus in front of A yy ). 
Disregarding the logarithmic terms, we find that the function A yy (T) is roughly a linear 
function of temperature for each type of pairing discussed. Furthermore, the interior part 
of the intervortex space ( T < pjv s (r) ), where one can use Eq. (|2lf) , results in K yy , which 
rises with magnetic field oc -R— 



In 



B C 2 



. These two regions compete with each other as 



it was above for k zz . But even in the case T ^ pfV s (r) large distances from vortex cores 



^ < r ^ t;^jB C 2/B dominate here leading to additional logarithmic factor, as compared to 
Eqs.([2l|), (|22|) . Possibly, this logarithmic factor is responsible for exceeding in several times 
of K yy as compared to n zz , which is experimentally observed in UPt^ 0,0. For quantitative 
consideration of this problem, one should evidently go beyond the approximation of locally 
quasiuniform superflow in describing the mixed state. 
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V. SCALING OF THE THERMAL CONDUCTIVITY AT LOW FIELDS AND 

LOW TEMPERATURES 



Thermodynamic and transport characteristics of superconductors with nodes may exhibit 
scaling behavior at low fields and low temperatures. This possibility for high-temperature 
superconductors was studied theoretically in |7|j2^-p9|. In particular, it was shown, that 
the electronic thermal conductivity of a two-dimensional <i-wave superconductor with four 
lines of nodes on a cylindrical Fermi surface may be represented as Kij ~ TFij(aT/B 1 ^ 2 ) 
| |28| |. One of the important assumptions underlying this result is that contributions to the 
thermal conductivity from quasiparticles in the intervortex space far enough from vortex 
cores dominate. In contrast to the thermodynamic quantities, for transport phenomena 
scattering processes are of importance even in clean limit. We note, that properties of 



disorder impurity potentials considered in Ref. [|8[] may be suitable for the Born scatterers, 
not in the unitary limit. 

Superconducting U Pt 3 essentially differs from the superconductors just mentioned above. 
It is three-dimensional hexagonal superconductor with the order parameter, which is believed 
to have both line of nodes on the equator of the Fermi surface and point nodes on its poles. 
While the strength of impurity scattering for high-temperature superconductors is not yet 
definitely determined, for heavy-fermion superconductors (in particular, for UPt 3 ) there are 
various experimental results and the physical background |f3ll,[40| indicating to the impurity 
scattering very close to the unitarity limit. Thus, the problem of scaling behavior of the 
thermal conductivity in superconductors like UPt 3 still has not been theoretically studied 
properly. 

At the same time the behavior quite close to the scaling one has been experimentally 
established recently for the thermal conductivity at low fields and low temperatures for 
superconducting UPt 3 both for the component n zz under the magnetic field parallel to z- 
axis and for n yy in the case of the field applied along ?/-axis @. These are just the cases 
which we consider in the article. Two forms of scaling are of interest for discussion of the 
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experimental data of Ref. [0: 
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(31) 
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(32) 



Kn(T,B = Q) 

In order to check whether our results are in agreement with any of these scaling behaviors, 
let us specify firstly power law dependences of the thermal conductivity upon temperature 
(under the condition 7 < T <C A max ) for three particular types of pairing in question in the 
absence of the magnetic field. 

According to Eq.(TJ), in order to estimate frequency dependence of the relaxation time 
in the unitarity limit, one should find the corresponding behavior of ImGo, ReGo- Contribu- 
tions to ImGo, ReGo from the line of nodes in the equatorial plane (subscript I), linear and 
quadratic point nodes at the poles of the Fermi surface (subscripts pi and p2 respectively) 
are as follows 

ReG ,/(^) : 
ReG , p i(^) 
ReG , p2 M 

One should emphasize that for all three cases we consider, the behavior of ReGo (i n contrast 
to ImGo) is governed, even at low frequencies, by the behavior of the order parameter not 
only near the nodes, but over the whole Fermi surface. So, the given expressions for ReGo 
are approximate ones and constants A u A p i, A p2 can't be determined unambiguously unless 
the behavior of the order parameter all over the Fermi surface is known. For the polar state 
A = Ai cos 6 and a spherical Fermi surface we get A\ = 1 from the comparison of flTT|) and 
(|33|) at u <C Ai- Let the order parameter in more general case have the form |A| = A Q f(x), 
where x = cos6> and f(x) m f'(0)x for \x\ < 1, \f(x)\ ~ \f(l)\(l^fx) for (l=Fx) < 1. Then 
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one can show for sufficiently small frequencies 
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where 



c 



1 1 

+ 



Chi A — lim ( [ 



i-g ds ln(2a|/'(0)|) hu>|f(l)|) \ 

lf(o)i 1 |/'(i)|' /(*) l/'(o)l J' {> 

So, in calculating A one has to integrate 1/ f(x) over the whole Fermi surface. 

On the other side, if one does not fix from the very beginning the behavior of the order 
parameter all over the Fermi surface but only near the nodes, constants of this origin may be 
considered as fitting parameters in comparison of the theoretical results with experimental 
data under the corresponding conditions. These fitting parameters, generally speaking, may 
manifest weak temperature dependence in the low temperature region, associated with the 
respective low-temperature dependence of the order parameter. Any particular choice of the 
parameters corresponds yet with a large number of basis functions of the representation, 
rather than only with a unique particular one. 

One can easily see from Eqs.(|33"D, fll4|) , that if only one kind of nodes is present, then at 
low frequencies 
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Low temperature dependence of the relaxation time of the form Tln 2 (A/T), found for the 
polar state in [|3l], is in agreement with the expression for t s j. 

For (1, i)-phases of E\ g - and i?2«-types of pairing with hybrid gap functions, having both 
the line of nodes and point nodes, the low-frequency dependences of the relaxation times are 
similar to the case of the polar state. Although the particular form of the order parameter 
all over the Fermi surface, for instance, point nodes result in the change of constants in the 
approximate expression for t s i, which then differ for each type of pairing. So, we write 
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where A e/ = 2AiA 2 /(Ai + 2A 2 ). 

Contributions from the line and point nodes to low- frequency behaviors of are 
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Substituting Eqs.( 
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we obtain the following leading low-temperature 



terms for k zz under the condition 7^T«T C 
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The first three results in Eq.([39|) (as well as in Eq . (|40j) ) are in agreement with those ob 



tained in [^TJ , where numerical coefficients were not specified. They concern the cases when 
only one kind of nodes is presented (both the line of nodes or one kind of point nodes, but 
not a hybrid gap function). The results in Eqs.(p9|), fl40| ) concerning hybrid gap functions, 
are new. The particular example of K 2Zj i 9 -component of the thermal conductivity demon- 
strates, that in the unitary limit in the presence of several kinds of nodes (that is for a hybrid 
gap function) the index of power law behavior of the low-temperature thermal conductivity 
may be greater than the least index among those taking place for superconductors with one 
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separate kind of the nodes ( in particular, in discussing of the E^-case - for superconductors 
with the linear point nodes ). 

Since the same relaxation time enters the expressions for both k zz and K yy , then for each 
particular type of pairing there are relations between the coefficients in those expressions, 
which govern the behavior of the anisotropy ratio K zz /K yy . For instance, as for E 2u type 
of pairing the both quantities I zz ,2u and I yVt 2u are proportional to u, the anisotropy ratio 
of the thermal conductivity in leading approximation is determined simply by the ratio 
Izz,2u/I yy ,2u = Ai(T)/A 2 (T), that is only weakly depends upon temperature. 

For E\ g case we have I Z z,ig oc u 2 , while Iyy t i g oc u. As a consequence, the anisotropy ratio 
of the thermal conductivity for E lg type of pairing essentially depends upon temperature: 



conductivity in UPt^ at low temperatures 7 ^ T C r c seems to permit one to distinguish 
between E 2u and E\ g representations in favor of (1, i)-phase of .E^-type of pairing. As it 
follows from Jl| for their particular clean samples, the temperature interval for which the low- 
temperature power law behavior of the thermal conductivity takes place is approximately 
0.07 < T/T c _ < 0.15. According to Eq. (|4~1|) , the anisotropy ratio should increase in more 
than 1.78 times, when the temperature changes from T = 0.07T C _ to T = 0.14T C _. This 
seems to be in a contradiction with the experiment (see, Fig. 2 in which shows the 
increase of the anisotropy ratio only on 10 — 12 per cent with the temperature change 
discussed. In other words, the fitting parameters available for the E\ a case in Eqs. (0), (|4"0"D 
allow to describe properly experimental results for the given temperature interval both for 
k zz or for Kyy, and not for both components simultaneously. 

While Eig type of pairing leads to noticeable overestimation of the low-temperature 




(41) 



Essential temperature dependence of the ratio for E\ g model was noticed earlier on the basis 
of numerical results in ]T8"| , f4"T| | . 

The analysis of experimental data of Ref. [|I|] on the anisotropy ratio of the thermal 
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dependence of the anisotropy ratio K zz /K yy in UPt 3 , within the framework of E 2u type of 
pairing we get, for the first sight, the underestimation of that dependence. Although for 
superconductors with nodes low-temperature deviation of the order parameter from its zero- 
temperature value is not exponentially small (as for s-wave isotropic case), but manifests 
power-law dependences, the temperature dependence of the quantity A 1 (T)/A 2 (T) is yet 
too weak in order to explain the observed change of the ratio on 10 — 12 per cent. However, 
the temperature dependence discussed can be described if one keeps in the expression for the 
anisotropy ratio, obtained within E 2u type of pairing, the first low-temperature correction 
to the leading term. For this purpose one should specify next terms in the expansion of the 
momentum dependence of the order parameter near nodes. We let |A| pa Ai\{-n/2 — 9) + 
L(n/2-6) 3 \ for \6 - (vr/2)| < 1, and |A| « A 2 (6 2 - D6 4 ) for \0\ < 1. Then the calculations 
of the anisotropy ratio under the condition 7 < T < T c result in 



is close to unity. 

We see, that this relation is sufficiently sensitive to the value of the coefficient D, while the 
effect of L is beyond the first correction to the leading term. Taking D as a fitting parameter, 
one can easily describe the observed change of the anisotropy ratio for 0.07 < T/T c _ < 0.15. 
For example, let A2 ~ 3Ai pa 12T C _, then we find D pa 5.2. Particular values for coefficients 
are obtained here for a spherical Fermi surface and they, of course, depend upon the form 
of the Fermi surface. It is of importance to emphasize, however, that the anisotropy of the 
Fermi surface doesn't change our qualitative conclusion itself. It is based on the fact, that 
K zz,i 9 l K yy,ig is roughly proportional to the temperature, while K ZZy2u / K yy,2u weakly depends 
on T, being at the same time sufficiently sensitive to the coefficient D in the next term of 
the expansion of momentum dependence of the order parameter near point nodes. Note 
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where the quantity 



S(T) 



ln 2 (A e/ /(5.6L4 2M T)) + 2.64 
ln 2 (A e/ /(4.67A 2u T)) + 2.67 



(43) 



25 



relatively large values of D, to which our analysis results in. The necessity for such a value 
of D should be taken into account when the particular basis function is chosen for the 
numerical study There is, of course, the problem of describing the temperature dependence 
of the anisotropy ratio at ultra low temperatures T < 7, which should take account of 
steeper slope of the curve, which the experiments show for K zz /K yy in the crossover between 
two regimes. This problem is not considered here. 

For higher temperatures the behavior of the anisotropy ratio becomes much more sen- 
sitive to the particular form of the order parameter all over the Fermi surface, not to its 
behavior mostly near nodes. As a consequence, the problem of discrimination between 
types of pairing in considering the anisotropy ratio of the thermal conductivity is essentially 



ambiguous in that case [18 



Furthermore, it follows from Eqs. (|2~l|) , (p2j), fl26|) , (|39|), ( [TOD ( in particular, from temper- 
ature dependences of A ZZ (T), A yy (T) ), that for temperatures in question scaling behavior 
(pl|) ( as well as ( |32|) ) is valid (in disregarding logarithmic terms) for both components 



k u (T, B) (i = x, y) and for both limiting cases 7 < T < T C ^B/B c2 , T C ^B/B c2 < T < T c 
just for ^^-representation. For E\ g type of pairing deviations from scaling are more notice- 



able for /^-component both in the temperature dominating region T C J B / B c2 <C T <C T c 



(if one uses the form ([31])), or in the field dominating region 7 < T < T c y B / B c2 (for 
the form (|32|) ). In order to distinguish between two models, considering the magnetic 
field dependence of k, one could introduce, along with the zero-field anisotropy ratio dis- 
cussed just now, the other characteristic as well. It is the anisotropy ratio of the form 
(k„(T, B) - k zz (T, B = 0))/ (K yy (T, B) - K yy (T, B = 0)) = A zz {T)/Ayy{T) under the tem- 
perature dominating condition. Analogously to the anisotropy ratio in zero field, quantity 



A ZZ {T) I 'A yy (T) , in accordance with Eqs.(|25|), (p9|), (pB]), is only weakly temperature depen- 
dent function within the framework of E 2u model, while it is roughly proportional to T for 
the Ei g type of pairing. The anisotropy of the Fermi surface doesn't change the conclusion. 
Unfortunately, for U Pt$ under the field B > B c \ low temperatures discussed do not belong 
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to the temperature dominating region T c y B / B c2 <T <T C . 

The other test for the discrimination between possible types of pairing in superconduct- 
ing UPts, which is based on the presence of superflow, is as follows. Since the only important 
qualitative difference between order parameters of E 2u and E\ g types of pairing is the multi- 
plicity of point nodes on the poles of the Fermi surface, it looks reasonable, for maximizing 
the difference between the effects, to consider k zz component of thermal conductivity under 
the presence of the supercurrent along the z-axis. In order not to mix various directions 
of the superflow in this case, let it be the uniform transport supercurrent in the absence of 
the magnetic field. Such a problem is associated with possible experiments on thin films 



or whiskers [42]. One can use for both types of pairing Eq.(15) for the relaxation time, 
inserting there 8q = 0. Further, under the given conditions 
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Inserting these results into (Q), we obtain under the temperature dominating conditions 
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and for the superflow dominating regime T <C PfV s <C A r 
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We see for the given direction of the supercurrent, that k zz is the nonmonotonous function 
upon the value of the supercurrent only in the case of E 2u type of pairing, while for the 
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E\ g model the thermal conductivity at low temperatures monotonously changes with the 
superflow velocity (see Fig. 3). This kind of experiments seems to be quite useful for the 
discrimination between candidates for the type of pairing in superconducting UPt%. 

One of the characteristic features of clean superconductors with nodes is the strong low 
energy dependence of the relaxation time, both in the Born approximation and in the unitary 
limit 19,0.11.13 ( see as wel1 Eqs.(D, flUD, (HI). (13) in the above )■ This is essential for 



finding the low temperature behavior of the thermal conductivity in clean limit. It is worth 
noting that in the unitary limit the relaxation time may manifest weak dependence upon 
energy for u/A max > 0.1, so that for not very clean samples with 0.lA max < 7 a model 



with energy independent relaxation time may present a reasonable approximation [f40| , |T9 
For these values of 7 the conditions 7 ^ T < T c , we are interested in in the article, are 
not satisfied. However, the low temperature region 7 < T < T c does exist for samples of 
Ref. |],|2[], for which their authors give the estimate 7 ~ 0.017-fT (one gets for this value 
7 « 0.038T C _ ~ 0.01A ma ,). 

We note, that in the presence of scattering by impurities, which is close to the unitarity 
limit, both the relaxation time and the electronic thermal conductivity of two-dimensional 
sufficiently clean d-wave superconductors would manifest analogous nonmonotonous depen- 
dences upon the magnetic field, which are found above for the three-dimensional supercon- 
ductors with nodes. 



VI. THERMAL CONDUCTIVITY DOMINATED BY IMPURITY-INDUCED 

BOUND STATES 

Impurity-induced bound states dominate the thermal conductivity under the condition 
T,VfV s 7- It is essential, that in this limit results turn out to be independent of the 
relation between T and VfV s . The spatial averaging over the intervortex space would con- 
tain both the region where VfV s < 7 and a part of the intervortex space where the opposite 
condition is valid ( since u/v s (£) ~ T c ). Hence impurity- induced bound states could influ- 
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ence essentially the spatially averaged thermal conductivity only for large enough values of 
7. So, for clean superconductors with sufficiently small 7 one can disregard the influence 
of impurity bound states on the magnetic field dependence of the thermal conductivity. 
Nevertheless, it is of interest to consider how the uniform superfluid flow influences the 
thermal conductivity dominated by impurity induced bound states. Below we show, that in 
the unitarity limit the thermal conductivity dominated by impurity induced bound states 
diminishes with increasing superflow field under various conditions. Hence, there is no non- 
monotonies in the superflow field dependence of the thermal conductivity at least unless the 
relation T,VfV s ^ 7 is broken. 

Excitation energy renormalized by impurities satisfies the relation ui(— u) = —u)*{u) 
T6| , |43[| . Thus, retaining for sufficiently small u a pair of terms in the expansion of uj{lo) in 



powers of a;, we have: 

Q(J) ~ 27 + auj , (48) 

where 7 and a are real positive parameters. Parameter a is a function of 7 and, in particular, 
in the unitary limit in the absence of the condensate flow field it takes the form 

-1/2 \ 



7 2 + |A| 2 



( 7 2 + 2|A| 2 ) ( 7 2 + |A| 



Sf 

For sufficiently small 7 one gets a = 1/2. In the presence of the superflow field parameter 7 
is a function of v s , and up to the first correction to its zero field value we get 

7 ~ 7 + bv 2 s . (50) 

To estimate b we make use of the results of Ref . || , and obtain that in the unitary limit 
this coefficient is negative and |fe| ~ f 2 /7o, at least within the logarithmic accuracy. 

We use the expansion of the integrand in Eq.(|I]) over a parameter (auj — v f\ s)^ / [ r y 2 + 
I A| 2 ), which is supposed to be small. Renormalized order parameter can be taken in this limit 
for u — 0, since we are interested in most important terms which are linear in temperature. 

29 



Retaining the corrections to the zero field term up to second order in the superflow field, we 
obtain 



TT NfT 
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v ( 7o 2 + |A| 2 ) 3 / 2 ( 7o 2 + |A| 2 ) 5 / 2 , 
In the presence of nodes of the order parameter only narrow regions on the Fermi surface 
near those nodes are of importance in Eq.flBTD under the condition 7 <C A max . This permits 
to consider contributions from each kind of nodes separately. Due to the same reason 
corrections containing v s and coming from the corresponding dependence of |A| in the first 
term in Eq . flBlf ) , may be disregarded in this case, since they would have an additional small 
factor ~ ( 7 o/A maa; ) 2 associated with small characteristic angular regions near the nodes. 

Let the order parameter have the line of nodes on the equator of a spherical Fermi surface: 
|A| = Ax 1 1 — 6\. Then we get from Eq. (|5lD the contribution from this line to the thermal 
conductivity K yy in the presence of the superflow, corresponding to the magnetic field along 
y- axis: 
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67o , 



(52) 



Contributions from point nodes on the poles ( both linear and quadratic ) are negligibly 
small for K yy . It is not the case for k zz . In the presence of magnetic field along z- axis we 
find from Eq. (|5l|) (within the logarithmic approximation) the following contribution from 
the line of nodes for this component of the thermal conductivity: 
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(53) 



At the same time the contribution from the linear point node is 



7T 



7o_ 

A 2 
^2 



A 2J 



(54) 



Making use of the estimation made above for the coefficient b, we have disregarded here the 
term, which is in ( 7 q/A2) 2 times less than the last term in ([54]). 
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Furthermore, for quadratic point node one gets 



K 22 



TT 2 N f V 2 f T 



2 2 i 
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(55) 

coincide with the ones 



. 2 24 7o A 2j 

In the absence of the superfluid flow field results (|5|) 
obtained in [pl|[n]]. According to Eqs. (|52"D and (|55|) , thermal conductivity diminishes with 
increasing superflow under the corresponding conditions. Negative value of b ensures the 
same qualitative conclusion for Eqs. (|53|) , (|54|) . Scaling relation for the thermal conductivity 
discussed in previous section is obviously broken in the limit T, v /v s < 7 under the conditions 
we consider. We note nonmonotonous superflow dependence of the thermal conductivity 
upon condensate flow field under the condition T < 7, since for VfV s < 7 the thermal 
conductivity diminishes while for VfV s ^> 7 it increases with increasing superflow field. 



VII. CONCLUSION 

In conclusion we have examined the possibility for nonmonotonous magnetic field de- 
pendence of the electronic thermal conduction along the magnetic field at low temperatures 
for sufficiently clean superconductors with nodes of the order parameter on the Fermi sur- 
face. We found that the contribution from low energy quasiparticles in the intervortex space 
is quite important in this respect and specific for superconductors with nodes. The effect 
comes from the influence of condensate flow field in the intervortex space on the scattering 
by nonmagnetic impurities of low-energy extended quasiparticles with momentum directions 
near nodes of the order parameter. The scattering is considered to be sufficiently close to 
the unitarity limit. We showed that the relaxation time at low energy is a nonmonotonous 
function upon the condensate flow field. 

Nonmonotonous magnetic field dependence of the thermal conductivity due to this con- 
tribution may take place for type II superconductors with large Ginzburg-Landau parameter 
within the temperature interval [3T c <J~^f^ < T < T c , where numerical factor j3 may be of 
the order of unity. These results are in a good qualitative agreement with recent experi- 
ments JHIl on the heavy- fermion superconductor UPt^. We explain the anomalously strong 
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influence of the magnetic field on the quasiparticle scattering, observing in this compound 
f4"2|| . We obtained as well that scaling behavior of the thermal conductivity with a single 
parameter well as weak low-temperature dependence of the anisotropy ratio 

Kzz/ftyy in zero field, are valid with logarithmic accuracy within the temperature interval 
7 ~ T <C T c for (l,z)-phase of .S^-type of pairing. Qualitatively it is quite close to that 
observed for UPt% in Ref. HHI. Under the same conditions E\ g model results in more no- 
ticeable deviations from the scaling, and in essential temperature dependence of the ratio 
k zz (T,B = Q)/nyy(T,B = 0). New test is proposed for discrimination between candidates 
for the type of pairing in UPts, based on the dependence of k zz upon the value of transport 
supercurrent flowing in thin films or whyskers along the hexagonal axis. 
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FIGURES 

FIG. 1. The relaxation time t s as a function of superffuid velocity v s , directed parallel to the 

basal plane, in the unitary limit for various energies uj/Aq(T) = 0.01 (1), 0.05 (2), 0.2 (3), 

1 (4), 5 (5) and three types of pairing: (a) - polar phase, A = Ao(T)cos# , (b) - (l,i)-state of 

£i 9 -pairing, |A| = Ao (T) | cos 9 | sin 9 , (c) - (l,i)-state of E^n-pairing, |A| = Ao (T) | cos 9 | sin 2 9 
FIG. 2. The thermal conductivity function of superfluid velocity v s , directed parallel 

to the basal plane, for various temperatures and three types of pairing: (a) - polar phase, (b) - 

£i 9 -pairing, (c) - -Ebu-pairing. The same basis functions are used as for Fig. 1 

FIG. 3. Normalized thermal conductivity k zz (T, v s )/k zz (T, 0) as a function of superfluid ve- 
locity v s , directed along z-axis, for unitary scatterers, various temperatures T/Aq(T) = 0.01 
(1), 0.025 (2), 0.05 (3), 0.1 (4), 0.2 (5) and two types of pairing: (a) - ^-pairing - 
|A| = A o (T)|cos0|sin0, (b) - £ 2u -pafring - |A| = A (T)| cos 9\ sin 2 9. 
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